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Abstract 

In this paper I developed a classical model of elementary particle that 
is associated with a membrane of finite size, surrounded by non-linear 
electromagnetic field. The form of local interaction which lead to bounded 
states of finite masses, charges and spins was constructed. To do this I 
added Kaluza-Klein Lagrangian on the surface, which is associated with 
extended particle, to quadratic in field potentials term (a la tachyon mass). 
This ensures that Lagrangian remains an even function of the field, but 
spontaneous symmetry breaking leads to nontrivial soliton-like solutions. 
I assumed that the particle has axial symmetry and that the surface has 
only one degree of freedom, i.e. is a disk with the radius determined from 
the equations of motion. The solution of system of two non-linear partial 
differential equations for the field potentials was obtained numerically by 
different methods. Several solutions with increasing orders of leading field 
multipoles and disk radius were obtained, and masses, electric charges and 
spins were calculated. In the framework of this model the ratio of a charge 
square to a double spin, i.e. the fine structure constant, which do not 
depend on parameters of the model, was calculated. 

PACS numbers: ll.10.Lm, 11.27.d, 12.90+b 

1 Introduction 

Physicists handle the laws of nature in two ways: they either explain them us- 
ing more fundamental hypotheses or just postulate them. The values of electric 
charges and spins of elementary particles are still almost postulated in the frame- 
work of quantum relativistic field theory, whereas their masses became a subject 
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of non-trivial calculations a long time ago. Wineberg-Salam model predicted not 
only existence and characteristics of intermediate bosons but allowed to calculate 
electroweak processes in tree-like approach. However, the calculation of radia- 
tion corrections is limited by the absence of experimental data on Higgs bosons' 
masses. Also further unifications of interactions such as supersymmetry [I] and 
great unification || were not supported by clear experimental observations so far, 
and electrodynamics are still leading in calculation precision of stationary states. 
Nevertheless electrodynamics can not be considered as an ideal model for the lep- 
ton world because of two problems: divergences in both classical and quantum 
theories, and the absence of calculations of some dimensionless quantities such 
as fine structure constant and lepton mass ratios. 

Non-linear extensions of electrodynamics were first proposed by Kaluza H], 
Klein |3|] as well as by Born and Infeld ||. Non-linearityleads to non-singular so- 
lution in case of point-like particle || , but the particle mass, charge and magnetic 
moment remain arbitrary. To change this situation one can introduce a funda- 
mental length into the model while preserving the local character of interaction, 
and construct a non-singular solution with finite size field sources. 

In the present paper I propose an extended particle model, which consists 
of only one real vector electromagnetic field. The model is based on bilinear 
Maxwell Lagrangian 

L M = ±- f(E 2 -B 2 )dV, (1) 

o 71 J 

and quadrolinear Kaluza-Klein Lagrangian (self-interaction of vector field) 

L "-"-shfJ (E - B ^ dV - < 2 > 

where constant (3 has dimension of field strength. 

The key assumption of this work is proportionality of the four-current on some 
surface associated with the particle to the four-potential of the field. Therefore 
the gauge invariance is broken on the surface. In order to generate this depen- 
dence in the form of Lagrange equations I add quadratic in potentials term on 
the surface to the Lagrangian 

Lsurface = T~~7 / (P* ~ A 2 ) dS, (3) 

where I is a fundamental length and p and A are scalar and vector potentials. 
The Lagrangian remains an even function of the field, and to provoke spontaneous 
symmetry breaking the sign of the last term must be the same as the sign used 
in the theory of vector tachyon field. 

The solution of the problem was found by using approximation of infinite 
number of degrees of freedom by the finite number. It appears that non-linearity 
of Kaluza-Klein interaction on the one hand and the competition of surface and 
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volume energy on the other hand lead to a spectrum of non-trivial stationary 
states with the finite values of all observables. This result was obtained by as- 
suming an axial symmetry of the whole system, and that the surface is a disk 
with only one degree of freedom — disk radius. The description of the field by 
only two degrees of freedom has already led to a set of stationary solutions. Then 
more field degrees of freedom were added (up to 1000) so that solutions achieved 
an asymptotic behavior. 

An important feature of the presented classical model is that the charges, 
magnetic moments and spins are no longer parameters of the Lagrangian, but 
arise due to a spontaneous symmetry breaking. These observables are propor- 
tional to some powers of two fundamental parameters of the Lagrangian (3 and 
I. To compare predictions of the model with experiment one can construct sev- 
eral dimensionless parameters, which do not depend on (3 and /. Therefore, in 
the following text I set these two parameters to unity without loss of generality. 
For each state I calculate Lande coefficient and the ratio of a charge square to a 
double spin e 2 /(2s), i.e. the fine structure constant. I also extract the ratios of 
masses, spins and magnetic moments of different states. 



2 Lagrangian in ellipsoidal coordinates 

Owing to an axial symmetry and particular choice of calibration the problem was 
reduced to solution of two non-linear partial differential equations for electric field 
potential p and the only non-zero component of vector-potential A = A^, which 
both depend on two spacial coordinates. I introduced Cartesian coordinates so 
that z is perpendicular to the disk plane, and then did the following transforma- 
tion to the ortogonal ellipsoidal coordinates 

x = p cos (0) , (4) 

y = p sin ((f)) , (5) 

sin(7r/2M;(M;)) ' v ; 



'1 - v(v) 2 

p = R- 



sin (ir/2w(w)) ' 

where R is the disk radius. This transformation reduced an infinite space to a 
unit square (w,v) (note that none of the functions depends on the third angular 
coordinate 0). I have concentrated on the mirror symmetry case L(—v) = L(v), 
and therefore evaluated all functions in the interval < v < 1. The transforma- 
tion |]-|6| can be tuned by two monotonically increasing functions w(w) and v(v ) 
(0 < w,w,v < 1). The simplest case is w(w) = w, v(v) = v, which we explore 
in subsection |3~T| . The particular choice of functions w(w) and v(v ) should not 
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affect the values of observables, but could be used to achieve better convergence 
to the asymptotic solution. I have also introduced contravariant potentials 

p=pR,A = ARh 4n (7) 

where = y/g^/R is dimensionless Lame coefficient, and "field strength" 

,_, d _ d _ . . 

= ~~dw P ' = ~dv P ' 

B W = ^-A,B V = ~A. (9) 
dv dw 

The boundary conditions on the axis v — 1 are A = with finite £!„, and 
at infinity w = are p = 0, y4 = 0. The boundary conditions in the disk plane 
outside the disk itself v — are E v = and B w = 0. Now we can write Lagrangian 
using new variables as 

L = L p + La + Lk-k, L p = L pv + L ps , La = La v + La s , (10) 

where subscripts v and s denote volume and surface terms, and 

R r l r\ 



L pv = -- [ [ (F w El + F v E 2 v )dvdw, (11) 
2 Jo Jo 

L Av = ^ C [\G w B 2 w + G v B 2 )dvdw, (12) 
2 Jo Jo 

L K -k = 7TE t t K ( B ■ E ? dv dw > (13) 
2 R Jo Jo 

L ps = Y £s»jPdv, (14) 
L As = -^- [ 1 S A A 2 dv, (15) 



2 Jo 



where 



F w — h v h^/hw, F v — h w h^jh v , 

K = \/{h v h w h,/,), (B-E) = (B w E w + E v B v ) 

S p = h v h^, S A = hy/htf,. 

Orbital contribution to the angular momentum is zero due to the axial sym- 
metry = 0. The ^-component of the spin is determined by the integral 

s = R 2 f 1 [\D w ^hl + D v ^h 2 v )Ah^dvdw, (16) 
Jo Jo dp dp 

where 

D W = E W + B w {B ■ E) K/F w , D v = E v + B v {B ■ E) K/F v . 
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3 Numerical solution 



In order to solve the model numerically, I approximated the fields by the finite 
number of degrees of freedom either by expansion of both potentials in series or by 
discretization of the unit square. The third method combines both approaches. 
In all methods the problem was reduced to a solution of the system of many 
non-linear algebraic equations. The system was solved by the iterative Newton's 
tangent method for many variables. I supposed that the approximate asymptotic 
solution was found if solutions converged and the energy values did not change 
significantly with increasing number of variables. To choose initial approximation 
I have used the property of Kaluza-Klein interaction (generally speaking, of any 
Lagrangian with only quadratic and fourth-order dependence on the fields), that 

L A + L K _ K = 0,L p + L K _ K = 0, (17) 

on exact solutions. To ensure positive values of energy H (in the following I use 
H = —L) one requires Ha — H p — —Hk-k > 0. 



3.1 Basic functions method 

The simplest method, which allows to find qualitative solutions with only two 
field degrees of freedom, is basic functions method. In this method one expresses 
potentials as a finite series of normalized functions with unknown amplitudes, so 
that Lagrangian becomes a function of these amplitudes. Calculation of bilinear 
part of Lagrangian requires evaluation of relatively small number of integrals, 
whereas number of integrals required for quadrolinear terms is proportional to 
fourth power of basic functions number. Therefore the advantages of the method 
can be seen only for the small number of correctly chosen functions. 

The natural choice of basic functions is solutions of Maxwell equations (fixing 
the choice of coordinate system so that w(w) = w and v(v) = v) with the first 
few multipole moments 

Pi = — , (18) 
~ _ (-7T w + sin (tt wj) (1 - v 2 ) 

— 1 + COS [TT W) 

_ _ ((4 + 2 cos (ir w)) ir w - Q sin (vrw)) (3/2w 2 - 1/2) 

P3 — 7—. 7 s , (20) 

— 1 + COS (TTW) 



which were normalized at v — and w = 1 to be n/2. The potentials become 

p(v,w) = J2P2k-i(v,w) a 2 fc_i, (21) 
k=i 
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A(v,w) = Y, A 2k{v,w)a 2kl (22) 

k=l 

where a\ = e/R, a 2 = ^/R 2 3/2, 03,... are unknown multipole amplitudes, e - 
charge, and fi - magnetic moment. 

Integration of bilinear part of energy leads to (for the first six multipoles): 

H pv = -R/2 (tt/2) (a\ + 16/5 a\ + 4096/729 a 2 5 ), (23) 

H Av = R/2 (tt/2) (4/3 a\ + 256/63 a\ + 16384/2475 a 2 6 ) , (24) 

H ps = R 2 /2 (it/2) 2 (1/2 a\ + 1/2 a\ + 1/2 a\ - 1/2 a x a 3 - 1/9 a x a 5 - 13/36 a 3 a 5 ), 

(25) 

H As = -R 2 /2 (tt/2) 2 (1/4 al+11/24 a^+29/60 Og-1/3 a 2 a 4 -l/12 a 2 a 6 -l/3 a 4 a 6 ). 

(26) 

Integration of quadrolinear part of energy is evaluated analytically over v and 
numerically over w. An approximate values of first 3x3 = 9 integrals are listed 
below 

H K _ K = 1/(2 R) (2.13a 2 a 2 + 6.19a 2 a 2 a 4 + 15.4 a\ a 2 + 22.2 a 2 a 2 

+ 16.7 ai a 2 a 3 a 4 + 121 a 2 a 2 — 12.0 ai a\ a 3 + 15.0 a 3 a\ a 4 — 98.5 a 2 a 4 a 2 ). (27) 

To find a simplest solution with only one electric and one magnetic multipoles, 
we note that the bilinear parts of the energies H p and Ha must be positive, which 
is achieved if 

1.27 < #1,^2 < 3.40 
4.07 < #3,^4 < 5.64 

7.15 < R 5 ,R 6 < 8.72 (28) 

10.3 < R 7 ,R 8 < 11.8 

13.4 < #9,^10 < 15.0 

where the subscripts denote the number of corresponding multipole. I considered 
values of radii which satisfy each line of eq. ^ separately in order to obtain 
solutions with minimal energy. First pair of multipoles is centrosymmetric electric 
field and dipolar magnetic moment. Therefore the energy f/i_ 2 is a function of 



three variables: two multipole amplitudes ai, a 2 and the radius R. From eqs. 23 



27 one gets: 



H 1>2 = -R/2 (1.57a 2 -2.09 a 2 )+i? 2 /2 (1.23 a 2 -0.617a 2 )-2.13 a\ a\/(2 R). (29) 

The values of these unknowns are determined from solution of the system of 
three algebraic equations, which are obtained from requiring that partial deriva- 
tives of the energy in respect to these variables are equal to zero. To obtain the 
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next state I chose two other degrees of freedom, which are quadrupole electric 
and octupole magnetic moments, and then solve analogous equations for H 3j 4, 
etc. 

Then the obtained solutions were used as a starting point for solution of La- 
grange equations with 8 degrees of freedom in order to increase the precision. 
This procedure did not lead to a convergent process for the first pair of multi- 
poles but did improve calculations for all others. The possible reason is that the 
determinant of the matrix of linearized equations does not have constant sign 
and can be close to zero, which lead to poor convergence. The physical reason 
of this problem is the areas with very strong electrical field where equation for 
magnetic field becomes hyperbolic rather than elliptic. 

The results of calculation for the first three pairs of leading multipoles are 
listed in Table 1 {m — H). 



Leading 


Included 


R 


m 


e 


multipoles 


multipoles 








1,2 


1..2 


2.81 


3.57 


3.25 


1,2 


1..6 








3,4 


3..4 


5.03 


0.431 




3,4 


1..8 


5.49 


0.366 


1.14 


5,6 


5..6 


8.05 


0.337 




5,6 


1..10 


8.05 


0.923 


2.05 



Table 1 



The table shows that as expected from eqs. |28| the disk size increases with 
increasing number of leading multipoles. Further increase of number of multipoles 
must be accompanied by simultaneous addition of new basic functions with non- 
Maxwellian it;- dependence for lower multipoles, which would lead to significant 
increase of number of integrals and calculation time. 



3.2 Lattice approach 

An alternative method is to define the potentials on some lattice and to assume 
linear interpolation between lattice points. Due to local character of interactions, 
presence of first derivatives only and linear interpolation assumption the energy 
depends only on the values of potentials at the neighboring sites, so that each 
algebraic equation contains limited number of variables (not more than 18), which 
does not increase with increasing number of lattice points. Therefore calculation 
time increases slower than the cube of number of points. 

Linear interpolation and integration by trapezium rule works well only for 
slowly varying functions. To increase smoothness of potentials I used specific 
choice of function v(v ), which "stretched" oscillations of the potentials near the 
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disk axis. In particular I used v(v) = 3/2 v — 1/2 v 3 and w(w) = w. The results 
of calculation for different number of lattice points are listed in Table 2. 



Leading 


Lattice 


R 


m 


e 


multipoles 


size 








3,4 


10 x 10 


6.01 


1.40 


3.55 


3,4 


12 x 12 


5.98 


1.43 


3.68 


3,4 


15 x 15 


5.96 


1.47 


3.81 


3,4 


20 x 20 


5.96 


1.52 


3.91 


3,4 


ooxoo 


5.96 


1.60 


4.04 


5,6 


12 x 12 


9.88 


1.83 


3.92 


5,6 


14 x 14 


9.62 


1.58 


3.72 


5,6 


18 x 18 


9.38 


1.41 


3.69 


5,6 


24 x 24 


9.23 


1.32 


3.68 


5,6 


ooxoo 


9.05 


1.22 


3.68 


7,8 


12 x 12 


15.0 


3.61 


5.89 


7,8 


14 x 14 


14.2 


2.44 


4.40 


7,8 


18 x 18 


13.4 


1.70 


3.96 


7,8 


24 x 24 


12.8 


1.40 


3.90 


7,8 


OO X oo 


12.1 


1.20 


3.88 



Table 2 



I used solutions of the previous section as a starting point for the lattice 
method. At the end of each part of the table extrapolated values for the infinite 
lattice are listed. For leading multipoles 1,2 solutions was found for lattice 6x6 
only. For leading multipoles 3,4 stable solutions exist for a large range of number 
of degrees of freedom, and observables quickly achieve their asymptotic values. 
Convergence gets worse with increasing order of leading multipole because higher 
order multipoles have more oscillations on the disk surface. In particular poten- 
tials of 3,4 multipole have two extremums, 5,6 — 3 extremums and 7,8 — 4 
extremums. 

3.3 Combined method 

To improve precision of calculations and the consistency of two previous methods 
I divided space into two parts: for large distances from the disk (w = 0..3/4) I 
used basic functions method with 2 x 5 = 10 functions, and in the vicinity of the 
disk (w = 3/4.. 1) I used lattice approach. The results of combined method are 
listed in Table 3. 
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Leading 


Lattice 


R 


m 


e 


mil H"lTl/~llf3c: 
111U.1 ul IJVJICo 


OlZjC 








Q A 
6,4 


1 9 v Q 
1Z X 6 


K Q7 

o.y i 


1/11 

1.41 


o.Oo 


1 A 
6,4 


1 fi. v A 
10 X 4 


o.yo 


1 /I7 
1.4 i 


Q 77 


Q A 
6,4 


ZU X 


K QP 

o.yo 


1 c;i 

1.01 


1 88 
O.OO 


6,4 


Zo X I 


^ QP 

o.yo 


1.00 


o.yz 


1 A 
6,4 


oo x oo/ 4 


K Q7 

o.y ( 


1.0 ( 


1 QA 

o.y4 


3,4 


ooxoo 


5.96 


1.60 


4.04 


0,\J 


1 P. y 4 


y.oi 


1 ^7 

l.O i 


^ 7n 

O. i u 


K P 
0,0 


zu X 


Q IP 

y .60 


1 /IQ 

i.4y 


1 Q1 

o.y i 


K p 
0,0 


OS v 7 
Zo X / 


O 91 

y.zi 


l.oO 


q 77 
6. I 1 


K P 
0,0 


QO v 8 
oZ X o 


(1 1 c 

y.io 


1 97 
l.Z ( 


q fit; 
O.OO 


K P 
0,0 


/in v m 

4U X 1U 


O 1 Q 

y.io 


1.Z0 


o.Oo 


K P 
0,0 


oo x oo/4 


o r\p 
y.uo 


1 1 o 

i. iy 


O.OO 


5,6 


OO X oo 


9.05 


1.22 


3.68 


7 8 


1 fi x 4 


1 S 7 




4 1 ^ 


7,8 


20 x 5 


13.1 


1.55 


3.87 


7,8 


28 x 7 


12.6 


1.32 


3.92 


7,8 


32 x 8 


12.5 


1.27 


3.91 


7,8 


40 x 10 


12.4 


1.21 


3.90 


7,8 


oo x oo/4 


12.1 


1.10 


3.90 


7,8 


oo x oo 


12.1 


1.20 


3.88 



Table 3 



The last lines of each part of this table repeats the asymptotic results of 
the previous section. The results of lattice and combined methods are in good 
agreement with each other. However, combined method is 10 times faster and 
therefore gives more precise results. The results for v = 2 v — v 2 and 40 x 10 
lattice differ from Table 3 only by a few percent even for 7,8 multipole. 

4 Conclusions 

I proposed a new concept of extended particle that is associated with a membrane 
of finite size. Charges and currents on this surface are chosen to be proportional 
to the corresponding potentials of the non-linear Kaluza-Klein electromagnetic 
field, so that the gauge invariance is broken on this surface. Then I assumed the 
axial and mirror symmetry of the system and existence of the only one surface 
degree of freedom: radius of the disk. Number of field degrees of freedom was 
also chosen to be finite. Described set of assumptions allowed me to develop 
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a technique for iterative calculation of the first three stationary states of the 
model. Then I demonstrated that increasing number of degrees of freedom (up 
to 1000) leads to the asymptotic solutions of Lagrange equations. Several states 
with increasing order of leading electric and magnetic multipoles and disk radii 
were constructed numerically using this method. Each state can be potentially 
associated with an elementary particle. 

The key result of the paper is that the values of all observables are finite. The 
masses were found to decrease weakly with increasing order of leading multipoles. 
The charges, magnetic moments and spins were also calculated from the model 
rather than postulated. It was not obvious that charges and spins of different 
states are identical or that their ratios are rational numbers. Calculations have 
shown that they are equal for the first three states within 10-20%: 



Leading 
multipoles 


R 


m 


e 


s 


l/a 


fl/flB 


3,4 


6.0 


1.6 


4.0 


12 


1.4 


1.0 


5,6 


9.1 


1.2 


3.6 


11 


1.5 


1.3 


7,8 


12.1 


1.1 


3.9 


12 


1.5 


1.1 



Table 4 



More detailed and extensive calculations are needed to make more precise 
conclusion. From the other side, the question about equality of charges and 
spins can be resolved by discovering some hidden symmetry of the described or 
modified interaction. It is also possible to change non-linearity to Born-Infeld 
type without adding any new parameters to the model. 

Masses, charges, magnetic moments and spins are proportional to some pow- 
ers of the two fundamental constants of the model with dimensions of length 
and field strength. However, there are two observables which do not depend on 
these constants — Lande coefficient and the ratio of charge square to a 

double spin, which is equal to the fine structure constant if one assumes that 
the spin is %/2. Lande coefficient was found about 1, and the fine structure con- 
stant — more than 100 times larger than expected. Therefore, calculated values 
of these parameters do not allow to associate discovered stationary states with 
charged leptons. Another interpretations (quarks) would require development of 
a calculation technique for bound states probably with non-additive charges and 
spins. 

The presented method allows to modify calculations for different non-linearity 
and/or boundary conditions on the surface without adding new model parame- 
ters. This will preserve the predictive ability of the approach and the existence 
of finite values of the observables. 
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